Frequency-domain method for joint equalization and decoding of space-time block codes

ABSTRACT

Space-time block coding is combined with single-carrier, minimum-mean-square-error, frequency-domain equalization for wireless communication through a medium characterized by multi-paths. The transmitter encodes incoming symbols into two or more streams that are transmitted over a corresponding number of transmitting antennas. The encoding employs modulo arithmetic. Decoding in the receiver proceeds by converting received signals to frequency domain, linearly combining the signals to separate contribution of the signals from the two or more transmitting antennas equalizing the separated signals, converting the equalized signals to time domain, and applying the converted signals to a decision circuit.

RELATED APPLICATIONS

This application is related to Provisional Application 60/282,634, filed Apr. 9, 2001.

BACKGROUND OF THE INVENTION

This relates to space-time coding, and more particularly, to frequency domain equalization in a receiver that receives signals from a multi-path channel.

M. V. Clark, in “Adaptive Frequency-Domain Equalization and Diversity Combining for Broadband Wireless Communications,” IEEE Journal on Selected Areas in Communications, pp. 1385-1395, October 1998, has shown that single-carrier minimum-mean-square-error frequency-domain (SC MMSE FDE) is an attractive equalization scheme for broadband wireless channels that are characterized by their long impulse response memories. Primarily the advantage of this approach accrues from using the computationally efficient Fast Fourier Transform (FFT), which compares favorably with time-domain equalization, where complexity grows exponentially with channel memory or require very long FIR filters to achieve acceptable performance. Further more, H. Sari et al, in “Transmission Techniques for Digital Terrestrial TV Broadcasting,” IEEE Communications Magazine, pp 100-190, February 1995, have shown that SC MMSE FDE has two main advantages over Orthogonal Frequency Division Multiplexing (OFDM); namely, lower peak-to-average ration, and reduced sensitivity to carrier frequency errors.

Diversity transmission using space-time block coding, for example, like disclosed in U.S. Pat. No. 6,185,258, has been proposed for several wireless standards because of the many attractive features of such coding. The SC MMSE FDE was combined with receive diversity by G. Kadel, in “Diversity and Equalization in Frequency Domain—A Robust and Flexible Receiver Technology for Broadband Mobile Communications Systems,” VTC, pp. 894-898, May 1997. There has also been some recent work on combining the scheme disclosed in the U.S. Pat. No. 6,185,258 with OFDM, as reported by Liu et al in “Decoding and Equalization of Unknown Multipath Channels Based on Block Precoding and Transmit-Antenna Diversity,” Asilomar Conf On Signals, Systems, and Computers, pp. 1557-1561, 1999. However, it has not been realized that various communication advantages result from combining space-time block coding with SC MMSE FDE.

SUMMARY OF THE INVENTION

An advance in the art is achieved by combining space-time block coding in a transmitter, for example, of a mobile unit, with SC MMSE FDE in a responsive receiver, for example in a base station. Illustratively, in a two transmitting antenna arrangement, the transmitter handles two blocks of N incoming symbols x_(a)(n) and x_(b)(n) at a time, encodes the block of N symbols into two streams, and transmits the two streams over two separate antennas in two consecutive frames, k, and k+1. In frame k, one antenna transmits symbols x₁ ^(k)(n), while the other antenna transmits symbols x₂ ^(k)(n). Symbols x₁ ^(k)(n) can, for example, equal symbols x_(a)(n), while symbols x₂ ^(k)(n) can equal x_(b)(n). In frame k+1, the one antenna transmits symbols x₁ ^(k+1)(n), while the other antenna transmits symbols x₂ ^(k+2)(n). The relationships adopted for x₁ ^(k+1)(n) and x₂ ^(k+2)(n) are x₁ ^(k+1) ^((n)=−{overscore (x)}) ₂ ^(k)((−n)^(N)) and x₂ ^(k+1)(n)=−{overscore (x)}₁ ^(k)((−n)_(N)) where {overscore (x)}₂ ^(k) is the complex conjugate of {overscore (x)}₂ ^(k). At the receiver, the received signal is transformed into the frequency domain, processed, equalized, transformed to time domain, and applied to a slicer to recover the transmitted sequences x₁ ^(k)(n) and x₂ ^(k)(n).

More specifically, the processing within the receiver multiplies Y, (which is the received signal that is transformed into the frequency domain) by Λ*. The matrix Λ* is the complex conjugate of Λ, where

${\Lambda = \begin{bmatrix} \Lambda_{1} & \Lambda_{2} \\ \Lambda_{2}^{*} & {- \Lambda_{1}^{*}} \end{bmatrix}},$ Λ₁ is a diagonal matrix with elements (q,q) that are equal to the q^(th) DFT coefficient of the transfer function between the one transmitting antenna and the receiver, and Λ₂ is a diagonal matrix with elements (q,q) that are equal to the q^(th) DFT coefficient of the transfer function between the other transmitting antenna and the receiver. Signal Y comprises signal Y^(k) for frame k, and signal Y^(k+1) for frame k+1. The frequency domain equalization is effected through multiplication of the signals forming vectors Y^(k) and Y^(k+1) by a diagonal matrix whose elements are related to the DFT coefficients between the one transmitting antenna and the receiver, and the other transmitting antenna and the receiver.

Extension of the principles disclosed herein to arrangements where several transmitting unit concurrently operate on the same carrier frequency is also disclosed.

BRIEF DESCRIPTION OF THE DRAWING

FIG. 1 shows a an illustrative arrangement of a two-antenna transmitter that is implemented in accordance with this disclosure, and a one antenna receiver that employs minimum-mean-square-error frequency-domain equalization (MMSE-FDE) equalization in accordance with this disclosure, with an interposed transmission channel having a memory of v symbol periods;

FIG. 2 shows the symbol stream that is generated by encoder 15 of FIG. 1; and

FIG. 3 presents an illustrative embodiment adapted for interference cancellation in an arrangement where several users employ the same carrier frequency to communicate with a receiver.

DETAILED DESCRIPTION

The following mathematical development focuses on a system having two transmit antennas and one receive antenna. It should be understood, however, that a skilled artisan could easily extend this mathematical development to more than two transmit antennas, and to more than one receive antenna.

FIG. 1 shows an arrangement that includes a transmitter with two transmit antennas 11 and 12, a receiver 20 with antenna 21, and a transmission channel therebetween that has a memory of v symbols periods. The transmission channel between antennas 11 and 21 can be represented as a finite impulse response filter with v+1 taps having values h₁(0), h₁(1), . . . h₁(v), forming vector h₁. Similarly the channel between antennas 12 and 21 can be represented by h₂. Within transmitter 10, encoder 15 outputs symbol blocks k=1, 2, 3, . . . , of two length N sequences; with one sequence, x₁ ^(k)(n), n=0, 1, 2, . . . N−1, being transmitted by antenna 11, and the other sequence, x₂ ^(k)(n), n=0, 1, 2, . . . N−1, being transmitted by antenna 12. In vector notation, the sequences can be expressed by x₁ ^(k) and x₂ ^(k), respectively. Because the transmission channel between transmitter 10 and receiver 20 has memory of v symbols, absent any protection, a receiver will commingle receptions from adjacent blocks. To avoid this, transmitter 10 includes a control signal that causes encoder 15 to precede each block with a string of v prefix symbols. In accordance with the principles disclosed herein, the symbols are selected to create a circulant transfer matrix between transmitter (for example, of a mobile unit) 10 and receiver to (for example, of a base station), as shown below.

It was discovered that advantageous results accrue from encoder 15 generating sequences such that the blocks of generated sequences are generated in pairs, k and k+1, where x ₁ ^(k+1)(n)=−{overscore (x)}₂ ^(k)((−n)_(N)) and x₂ ^(k+1)(n)={overscore (x)}₁ ^(k)((−n)_(N))  (1) for n=0, 1, . . . N−1 and k=0, 2, 4, . . . , where {overscore (x)}_(a) denotes complex conjugation of x_(a), (subscript a being either 1 or 2), and (.)_(N) denoted modulo-N operation. In other words, when encoder 15 generates the information sequence

{x(0), x(1), x(2), . . . x(N−2)x(N−1)},

for antenna 11 for frame k, it concurrently generates the information sequence

{{overscore (x)}(0), {overscore (x)}(N−1),{overscore (x)}(N−2), . . . {overscore (x)}(2),{overscore (x)}(1),}

for antenna 12 for or frame k+1. From equation (1) and the fact that X_(a) ^(k)=Qx_(a) ^(k), where Q is the orthogonal Discrete Fourier Transform (DFT) matrix whose p,q) element is given by

${{Q\left( {p,q} \right)} = {\frac{1}{\sqrt{N}}{\mathbb{e}}^{{- j}\frac{2\pi}{N}p\; q}}},$ and 0≦p,q≦N−1, is it follows that X ₁ ^(k+1)(m)=−{overscore (X)}₂ ^(k)(m) and X₂ ^(k+1)(m)={overscore (X)}₁ ^(k)(m)  (2)

As will become clear from the disclosure below, the relationship of equation (2) allows for frequency domain decoupling in a receiver the contributions of sequences x₁ ^(k)(n) and x₂ ^(k)(n). Accordingly, for purposes of this disclosure, such encoding—the encoding of equation (1) being one example thereof, is termed “frequency domain decouplable contributions” (FD-DC) encoding. It is encompassed by a class of orthogonal designs, as disclosed, for example, in U.S. Pat. No. 6,088,408, issued Jul. 11, 2000

FIG. 2 presents the image of a sequence that encoder 15 delivers to antenna 11, with the prefix string consisting symbols x_(i)(−v) through x₁(−1), followed by symbols x₁(0) through x₁(N−1). Ignoring the subscript, the signals received at antenna 21 in response to signals transmitted by antenna 11, starting with signal y(0), which corresponds to the received signal at the time the encoded signal x(0) is transmitted by antenna 11, can be expressed by the following matrix equation:

$\begin{matrix} {\begin{bmatrix} {y(0)} \\ {y(1)} \\ \vdots \\ {y(v)} \\ {y\left( {v + 1} \right)} \\ \vdots \\ {y\left( {N - 1} \right)} \end{bmatrix} = {\begin{bmatrix} {h(v)} & {h\left( {v - 1} \right)} & \cdots & {h(1)} & {h(0)} & 0 & \cdots & 0 \\ 0 & {h(v)} & \cdots & {h(2)} & {h(1)} & {h(0)} & \cdots & 0 \\ \vdots & \vdots & ⋰ & \vdots & \vdots & \vdots & ⋰ & \vdots \\ 0 & 0 & \cdots & 0 & {h(v)} & {h\left( {v - 1} \right)} & \cdots & 0 \\ 0 & 0 & \cdots & 0 & 0 & {h(v)} & \cdots & 0 \\ \vdots & \vdots & ⋰ & \vdots & \vdots & \vdots & ⋰ & \vdots \\ 0 & 0 & \cdots & 0 & 0 & 0 & \cdots & {h(0)} \end{bmatrix}\begin{bmatrix} {x\left( {- v} \right)} \\ {x\left( {{- v} + 1} \right)} \\ \vdots \\ {x\left( {- 1} \right)} \\ {x(0)} \\ {x(1)} \\ \vdots \\ {x\left( {N - 1} \right)} \end{bmatrix}}} & (3) \end{matrix}$ If the prefix sequence is chosen so that x(−i) x(N−i) for i=1, 2, . . . v, then equation (3) can be replaced by

$\begin{matrix} {\begin{bmatrix} {y(0)} \\ {y(1)} \\ \vdots \\ {y(v)} \\ {y\left( {v + 1} \right)} \\ \vdots \\ {y\left( {N - 1} \right)} \end{bmatrix} = {{\begin{bmatrix} {h(0)} & 0 & 0 & \cdots & {h(2)} & {h(1)} \\ {h(1)} & {h(0)} & 0 & \cdots & {h(3)} & {h(2)} \\ \vdots & \vdots & \vdots & ⋰ & \vdots & \vdots \\ {h\left( {v - 1} \right)} & {h\left( {v - 2} \right)} & {h\left( {v - 3} \right)} & \cdots & 0 & {h(v)} \\ {h(v)} & {h\left( {v - 1} \right)} & {h\left( {v - 2} \right)} & \cdots & 0 & 0 \\ \vdots & \vdots & \vdots & ⋰ & \vdots & \vdots \\ 0 & 0 & 0 & \cdots & {h(1)} & {h(0)} \end{bmatrix}\begin{bmatrix} {x(0)} \\ {x(1)} \\ {x(2)} \\ \vdots \\ {x\left( {N - 2} \right)} \\ {x\left( {N - 1} \right)} \end{bmatrix}}.}} & (4) \end{matrix}$ Returning to the use of subscripts to indicate the transmitting antenna, and superscripts to indicate symbol blocks, and recognizing that antenna 21 receives signals transmitted by antenna 12 as well as by antenna 11, we get y ^(k) =H ₁x₁ ^(k) +H ₂x₂ ^(k) +n ^(k)  (5) where

$\begin{matrix} {H_{a} = \begin{bmatrix} {h_{a}(0)} & 0 & \cdots & 0 & 0 & {h_{a}(v)} & \cdots & {h_{a}(1)} \\ {h_{a}(1)} & {h_{a}(0)} & \cdots & 0 & 0 & 0 & \cdots & {h_{a}(2)} \\ \vdots & \vdots & ⋰ & \vdots & \vdots & \vdots & ⋰ & \vdots \\ {h_{a}(v)} & {h_{a}\left( {v - 1} \right)} & \cdots & {h_{a}(0)} & 0 & 0 & \cdots & 0 \\ 0 & {h_{a}(v)} & \cdots & {h_{a}(1)} & {h_{a}(0)} & 0 & \cdots & 0 \\ 0 & 0 & \cdots & {h_{a}(2)} & {h_{a}(1)} & {h_{a}(0)} & \cdots & 0 \\ \vdots & \vdots & ⋰ & \vdots & \vdots & \vdots & ⋰ & \vdots \\ 0 & 0 & \cdots & 0 & {h_{a}(v)} & {h_{a}\left( {v - 1} \right)} & \cdots & {h_{a}(0)} \end{bmatrix}} & (6) \end{matrix}$ and the number of rows and columns in H_(a) is N.

A matrix with the structure of H_(a) is known as a circulant matrix, and it can be shown to have the eigen-decomposition H _(a) =Q ⁻¹Λ_(a) Q  (7) where (.)* denotes complex-conjugate transpose, Q⁻¹ is the inverse DFT matrix, and Λ_(a) is a diagonal matrix whose (q,q) element is equal to the q^(th) DFT coefficient of ha. It may be noted that Q⁻¹=Q* which, of course simplifies processing since it is not necessary to compute the inverse of the matrix, and it is assumed that the characteristics of the transmission medium between transmitter 10 and receiver 20 are knows by means of any one of knows processes. See, for example, U.S. patent application Ser. No. 09/956,648 filed Sep. 20, 2001 by the inventor of this application. Accordingly, it is presumed that the DFT coefficients of h_(a) are known and, therefore, Λ_(a) is known.

It may be noted also that the above assumes that the transmission channel between transmitter 10 and receiver 20 does not change during the course of computations that are relevant to this disclosure, which is the time of two consecutive blocks.

Applying the output signals of circuit 22 to FFT circuit 23 to transform the signals of equations (5) into the frequency domain yields

$\begin{matrix} \begin{matrix} {{Y^{k} \equiv {Q\; y^{k}}} = {{{Q\left( {Q^{- 1}\Lambda_{1}Q} \right)}x_{1}^{k}} + {{Q\left( {Q^{- 1}\Lambda_{2}Q} \right)}x_{2}^{k}} + {Q\; n^{k}}}} \\ {= {{\Lambda_{1}Q\; x_{1}^{k}} + {\Lambda_{2}Q\; x_{2}^{k}} + {Q\; n^{k}}}} \\ {= {{\Lambda_{1}X_{1}^{k}} + {\Lambda_{2}X_{2}^{k}} + {N^{k}.}}} \end{matrix} & (8) \end{matrix}$ and Y ^(k+1)=Λ₁ X ₁ ^(k+1)+Λ₂ X ₂ ^(k+1) +N _(k).  (9)

By taking advantage of the particular encoding by encoder 15, as expressed in equation (2), one can note that equation (9) can be rewritten and expressed together with equation (8) by

$\begin{matrix} {{{Y \equiv \begin{bmatrix} Y^{k} \\ {\overset{\_}{Y}}^{k + 1} \end{bmatrix}} = {{\begin{bmatrix} \Lambda_{1} & \Lambda_{2} \\ \Lambda_{2}^{*} & {- \Lambda_{1}^{*}} \end{bmatrix}\begin{bmatrix} X_{1}^{k} \\ X_{2}^{k} \end{bmatrix}} + \begin{bmatrix} N^{k} \\ {\overset{\_}{N}}^{k} \end{bmatrix}}}{Y \equiv {{\Lambda\; X} + {N.}}}} & (10) \end{matrix}$ Noting that Λ is an orthogonal matrix, the two signals X₁ ^(k) and X₂ ^(k) can be decoupled by simply multiplying both sides of equation (10) by Λ*, resulting in

$\begin{matrix} \begin{matrix} {{\overset{\sim}{Y} \equiv \begin{bmatrix} {\overset{\sim}{Y}}^{k} \\ {\overset{\sim}{Y}}^{k + 1} \end{bmatrix}} = {\Lambda^{*}Y}} \\ {= {{\begin{bmatrix} {{\Lambda_{1}\Lambda_{1}^{*}} + {\Lambda_{2}\Lambda_{2}^{*}}} & 0 \\ 0 & {{\Lambda_{1}\Lambda_{1}^{*}} + {\Lambda_{2}\Lambda_{2}^{*}}} \end{bmatrix}\begin{bmatrix} X_{1}^{k} \\ X_{2}^{k} \end{bmatrix}} + {\begin{bmatrix} {\overset{\sim}{N}}^{k} \\ {\overset{\sim}{N}}^{k + 1} \end{bmatrix}.}}} \end{matrix} & (11) \end{matrix}$ Expressed more succinctly,

$\begin{matrix} {{{\overset{\sim}{Y} \equiv \begin{bmatrix} {\overset{\sim}{Y}}^{k} \\ {\overset{\sim}{Y}}^{k + 1} \end{bmatrix}} = {{\begin{bmatrix} \overset{\sim}{\Lambda} & 0 \\ 0 & \overset{\sim}{\Lambda} \end{bmatrix}\begin{bmatrix} X_{1}^{k} \\ X_{2}^{k} \end{bmatrix}} + \begin{bmatrix} {\overset{\sim}{N}}^{k} \\ {\overset{\sim}{N}}^{k + 1} \end{bmatrix}}},} & (12) \end{matrix}$ from which we get

$\begin{matrix} {{{\overset{\sim}{Y}}^{k} = {{\overset{\sim}{\Lambda}X_{1}^{k}} + {\overset{\sim}{N}}^{k}}}{{\overset{\sim}{Y}}^{k + 1} = {{\overset{\sim}{\Lambda}X_{2}^{k}} + {\overset{\sim}{N}}^{k - 1}}}} & (13) \end{matrix}$ where {tilde over (Λ)}=(Λ₁Λ*+Λ₂Λ*) in an N×N diagonal matrix.

Signals {tilde over (Y)}^(k+1) and {tilde over (Y)}^(k) are computed in circuit 24 from

$\begin{matrix} {{\begin{bmatrix} {\overset{\sim}{Y}}^{k} \\ {\overset{\sim}{Y}}^{k + 1} \end{bmatrix} = {{\Lambda^{*}Y} = {\begin{bmatrix} \Lambda_{1}^{*} & \Lambda_{2} \\ \Lambda_{2}^{*} & {- \Lambda_{1}} \end{bmatrix}\begin{bmatrix} Y^{k} \\ {\overset{\_}{Y}}^{k + 1} \end{bmatrix}}}},} & (14) \end{matrix}$ which corresponds to {tilde over (Y)} ^(k)=Λ₁*Y^(k)+Λ₂{tilde over (Y)}^(k+1)  (15) and {tilde over (Y)} ^(k+1)=Λ₂*Y^(k)−Λ₁{tilde over (Y)}^(k+1),  (16) where signal Y^(k) is the output of FFT circuit 23 at block k, and signal {tilde over (Y)}^(k+1) is the complex conjugate of the output of FFT circuit 23 at block k+1.

In accord with the principles disclosed herein, for the type transmitter structure disclosed above, a single-carrier minimum-mean-square-error frequency-domain equalization (SC-MMSE-FDE) is obtained with an N×N diagonal matrix applied to the signals of equations (15) and (16), where the elements of this matrix are given by

$\begin{matrix} {{W\left( {i,i} \right)} = {\frac{1}{{\overset{\sim}{\Lambda}\left( {i,i} \right)} + \frac{1}{S\; N\; R}}.}} & (17) \end{matrix}$ where

${{S\; N\; R} \equiv \frac{\sigma_{x}^{2}}{\sigma_{n}^{2}}},$ , under the assumption that both the noise and the signal are “white.” It is noted that the noise vector n and its DFT N have the same variance, since the are related through the orthonormal transformation, Q.

Applying the output signals obtained by combiner 24 to SC-MMSE-FDE circuit 25, wherein the matrix of equation (17) pre-multiplies the signals of equations (15) and (16), results in

$\begin{matrix} \begin{matrix} {Z^{k} = {W{\overset{\sim}{Y}}^{k}}} \\ {= {{W\overset{\sim}{\Lambda}X_{1}^{k}} + {W{\overset{\sim}{N}}^{k}}}} \end{matrix} & (18) \end{matrix}$ where W{tilde over (Λ)} is a diagonal N×N matrix with elements

$\frac{\overset{\sim}{\Lambda}\left( {i,i} \right)}{{\overset{\sim}{\Lambda}\left( {i,i} \right)} + \frac{1}{S\; N\; R}},$ , and Z ^(k+1) =W{tilde over (Λ)}X ₁ ^(k+1) +WÑ ^(k+1)  (19) Applying equation (18) and equation (19) signals of circuit 25 to inverse FFT transform circuit 26 yields z ^(k) =Q ⁻¹ Z ^(k) =Q*X ^(k) =Q*W{tilde over (Λ)}X ₁ ^(k) +Q*WÑ ^(k) =Q*W{tilde over (Λ)}Qx ₁ ^(k) +Q*WQñ ^(k) ={tilde over (H)}x ₁ ^(k) +ñ  (20) and z ^(k+1) ={tilde over (H)}x ₂ ^(k) +ñ,  (21) where {tilde over (H)}=Q*W{tilde over (Λ)}Q has the same form as H in equation (7). Note that since the noise statistics are the same in frames k and k+1, equations (20) and (21) replace both Q*WQñ^(k) and Q*WQñ^(k+1) with ñ. Also note the {tilde over (H)} approaches the identity matrix as the value of SNR increases. For this reason,

Thus, the signals developed by circuit 26 correspond to the signals x₁ ^(k) and x₂ ^(k) that are filtered and equalized by matrix {tilde over (H)}, and augmented by noise. It is noted that the matrix {tilde over (H)} approaches the identity matrix as the value of SNR increases. For this reason a simple “slicer” detector can be used to detect the signals x₁ ^(k) and x₂ ^(k).

Accordingly, the signals developed by circuit 26 are applied to detector 27, which is a conventional slicer, yielding the signals x₁ ^(k) and x₂ ^(k), from which the sequence of transmitted symbols is obtained.

Expanding on the principles disclosed herein, FIG. 3 depicts an arrangement where M units (two shown) employ the same channel (i.e., the same carrier frequency) to communicate with receiver 20, which has M antennas (two shown). Corresponding to such an arrangement, equation (10) is rewritten as

$\begin{matrix} {{\begin{bmatrix} Y_{21} \\ Y_{22} \end{bmatrix} = {{\begin{bmatrix} \Lambda_{10 - 21} & \Lambda_{30 - 21} \\ \Lambda_{10 - 22} & \Lambda_{30 - 22} \end{bmatrix}\begin{bmatrix} X_{10} \\ X_{30} \end{bmatrix}} + \begin{bmatrix} N_{21} \\ N_{22} \end{bmatrix}}},} & (22) \end{matrix}$ where Y₂₁ represents the signals received at antenna 21, X₁₀ represents signals transmitted by unit 10, Λ₁₀₋₂₁ represents the channel between unit 10 and antenna 21, X₃₀ represents signals transmitted by unit 30, Λ₃₀₋₂₁ represents the channel between unit 30 and antenna 21, and N₂, is the noise received at antenna 21. Similarly, Y₂₂ represents the signals received at antenna 22, Λ₃₀₋₂₂ represents the channel between unit 30 and antenna 22, Λ₁₀₋₂₂ represents the channel between unit 10 and antenna 22, and N₂₂ is the noise received at antenna 22.

The signals of the two units can be decoupled in each frame by applying the following linear, zero-forcing interference cancellation:

$\begin{matrix} \begin{matrix} {\begin{bmatrix} Y_{10}^{k} \\ Y_{30}^{k} \end{bmatrix} = {\begin{bmatrix} I & {{- \Lambda_{30 - 21}}\Lambda_{10 - 22}^{- 1}} \\ {{- \Lambda_{30 - 22}}\Lambda_{10 - 21}^{- 1}} & I \end{bmatrix}\begin{bmatrix} Y_{21}^{k} \\ Y_{22}^{k} \end{bmatrix}}} \\ {= {{{\begin{bmatrix} I & {{- \Lambda_{30 - 21}}\Lambda_{30 - 22}^{- 1}} \\ {{- \Lambda_{10 - 22}}\Lambda_{10 - 21}^{- 1}} & I \end{bmatrix}\begin{bmatrix} \Lambda_{10 - 21} & \Lambda_{30 - 21} \\ \Lambda_{10 - 22} & \Lambda_{30 - 22} \end{bmatrix}}\begin{bmatrix} X_{10}^{k} \\ X_{30}^{k} \end{bmatrix}} + \begin{bmatrix} {\overset{\sim}{N}}_{21}^{k} \\ {\overset{\sim}{N}}_{22}^{k} \end{bmatrix}}} \\ {= {{\begin{bmatrix} {\overset{\sim}{\Lambda}}_{10 - 21} & 0 \\ 0 & {\overset{\sim}{\Lambda}}_{10 - 22} \end{bmatrix}\begin{bmatrix} X_{10}^{k} \\ X_{30}^{k} \end{bmatrix}} + \begin{bmatrix} {\overset{\sim}{N}}_{21}^{k} \\ {\overset{\sim}{N}}_{22}^{k} \end{bmatrix}}} \end{matrix} & (23) \end{matrix}$ where Λ₁₀₋₂₁=Λ₁₀₋₂₁−Λ₃₀₋₂₁Λ₃₀₋₂₂ ⁻¹Λ₁₀₋₂₂ and {tilde over (Λ)}₁₀₋₂₂=Λ₃₀₋₂₂−Λ₁₀₋₂₂Λ₁₀₋₂₁ ⁻¹Λ₃₀₋₂₁. The critical observation to make here is that both {tilde over (Λ)}₁₀₋₂₁ and {tilde over (Λ)}₁₀₋₂₂ are orthogonal matrices like matrix Λ in equation (10) and, therefore, decoding proceeds as described above for the single unit case, and the full diversity gain is guaranteed for both users. That is, combiner 28 in FIG. 3, which is a modified version of combiner 24 in FIG. 1, first performs the processing called for by equation 23, and then proceeds as described above. 

1. A receiver comprising: a time-domain to frequency-domain converter responsive to a signal received by an antenna in frames k and k+1, for developing signals Y^(k) in frame k and signals Y^(k+1) in frame k+1; a linear combiner for creating a first linear combination signal, {tilde over (Y)}^(k), from signals related to Y^(k) and Y^(k+1), and a second linear combination signal, {tilde over (Y)}^(k+1), from signals related to Y^(k) and Y^(k+1), where said first linear combination is different from said second linear combination; an equalizer that pre-multiplies signal {tilde over (Y)}^(k) by a diagonal matrix W to form signal {tilde over (X)}^(k), and pre-multiplies signal {tilde over (Y)}^(k+1) by said diagonal matrix W to form signal {tilde over (Z)}^(k+1); a frequency-domain to time-domain converter for converting signals {tilde over (Z)}^(k) and {tilde over (Z)}^(k+1) to time-domain signals; and a slicer responsive to said time domain signals.
 2. The receiver of claim 1 where said time-domain to frequency-domain converter implements a Fast Fourier Transform algorithm.
 3. The receiver of claim 1 where said frequency-domain to time-domain converter implements an inverse Fast Fourier Transform algorithm.
 4. The receiver of claim 1 where said linear combiner, in creating signal {tilde over (Y)}^(k) from component signals related to Y^(k) and Y^(k+1), multiplies at least one of said component signals by a diagonal matrix.
 5. The receiver of claim 1 where said linear combiner, in creating signal {tilde over (Y)}^(k) from component signals related to Y^(k) and Y^(k+1), multiplies each of said component signals by a different diagonal matrix.
 6. The receiver of claim 1 where said linear combiner, in creating signal {tilde over (Y)}^(k) from component signals related to Y^(k) and Y^(k+1), employs diagonal matrices Λ₁ and Λ₂ where diagonal matrix Λ₁ is related to characteristics of transmission medium between a first antenna of a transmitter of signals received by said receiver, and Λ₂ is related to characteristics of transmission medium between a first antenna of a transmitter of signals received by said receiver.
 7. The receiver of claim 6 where said linear combiner, in creating signal {tilde over (Y)}^(k+1) from component signals related to Y^(k) and Y^(k+1), employs diagonal matrices that are related to said matrices Λ₁ and Λ₂ through operations taken from a set that includes negations and complex conjugations.
 8. The receiver of claim 1 where said linear combiner creates signal {tilde over (Y)}^(l)=Λ₁ *Y ^(k)+Λ₂ {overscore (Y)} ^(k+1), and signal {tilde over (Y)}^(k+1)=Λ₂ *Y ^(k)−Λ₁ {overscore (Y)} ^(k+1), where {overscore (Y)}^(k+1) is a complex conjugate of Y^(k+1).
 9. The receiver of claim 8 where elements of said diagonal matrix W are related to matrices Λ_(1 and Λ) ₂.
 10. The receiver of claim 8 where said diagonal matrix W has elements ${{W\left( {i,i} \right)} = \frac{1}{{\overset{\sim}{\Lambda}\left( {i,i} \right)} + \frac{1}{S\; N\; R}}},{{{where}\mspace{14mu}{\overset{\sim}{\Lambda}\left( {i,i} \right)}} = {{{\Lambda_{1}\left( {i,i} \right)}{\Lambda_{1}^{*}\left( {i,i} \right)}} + {{\Lambda_{2}\left( {i,i} \right)}{\Lambda_{2}^{*}\left( {i,i} \right)}}}},$ and (.)* represents a complex conjugate operation, and SNR is a computed value.
 11. A receiver comprising: a time-domain to frequency-domain converter responsive to a signal received by an antenna in frames k, k+1, . . . k+m, where m is a selected integer greater than 0, for developing signals Y^(k), Y^(k+1), . . . Y^(k+m), in frames k, k+1, . . . . k+m, respectively; a linear combiner for creating signals {tilde over (Y)}^(k), {tilde over (Y)}^(k+1), . . . {tilde over (Y)}^(k+m) from linear combinations of signals related to Y^(k), Y^(k+1), . . . Y^(k+m); an equalizer that pre-multiplies each signal {tilde over (Y)}^(j),j=k, k+1, . . . k+m by a diagonal matrix W to form signals {overscore (Z)}^(l), j=k, k+1, . . . k+m; a frequency-domain to time-domain converter for converting signals {tilde over (Z)}^(j) to time domain signals; and a slicer responsive to said time domain signals.
 12. The receiver of claim 11 where said signals related to signals Y^(k), Y^(k+1), . . . y^(k+m) are related to said signals Y^(k), Y^(k+1), . . . Y^(k+m) through operations from a set that includes negations and complex conjugations.
 13. A receiver comprising: p antennas, where p is an integer greater than 1; a time-domain to frequency-domain converter responsive to a signal received by each of said antennas in frames k, k+1, . . . k+m, where m is a selected integer greater than 0, for developing signals Y_(j) ^(k), Y_(j) ^(k+1), . . . y_(h) ^(k+m), in frames k, k+1, . . . k+m, respectively, where subscript j identifies a j^(th) antennas of said p antennas; a linear combiner for creating groups of signals {tilde over (Y)}_(n) ^(k), {tilde over (Y)}_(n) ^(k+1), . . . {tilde over (Y)}_(n) ^(k+m) for each value of subscript j=1, 2, . . . p, from linear combinations of signals related to said signals {tilde over (Y)}_(n) ^(k), {tilde over (Y)}_(n) ^(k+1), . . . {tilde over (Y)}_(n) ^(k+m), when n is an index designating a transmitting unit that supplies signals to said p antennas; an equalizer that pre-multiplies each signal {tilde over (Y)}_(n) ^(q), q=k, k+1, . . . k+m by a diagonal matrix W to form signals {tilde over (Z)}_(n) ^(q), q=k, k+1, . . . k+m; a frequency-domain to time-domain converter for converting signals {tilde over (Z)}_(n) ^(q) to time-domain signals; and a slicer responsive to said time domain signals.
 14. The receiver of claim 13 where p=2, and where said linear combiner obtains signals {tilde over (Y)}_(n) ^(k) and {tilde over (Y)}_(n) ^(k+1) by computing $\begin{bmatrix} {\hat{Y}}_{1}^{k} \\ {\hat{Y}}_{2}^{k} \end{bmatrix} = {\begin{bmatrix} I & {{- \Lambda_{2 - 1}}\Lambda_{1 - 2}^{- 1}} \\ {{- \Lambda_{2 - 2}}\Lambda_{1 - 1}^{- 1}} & I \end{bmatrix}\begin{bmatrix} Y_{1}^{k} \\ Y_{2}^{k} \end{bmatrix}}$ where Ŷ₁ ^(k) represents signal received at said receiver, in frame k, from transmitting unit 1, and Ŷ₂ ^(k) represents signal received at said receiver, in frame k, from transmitting unit 1, Λ₁₋₁ is a diagonal matrix representing transmission medium between transmitting unit 1 and a first one of said two antennas, Λ₂₋₁ is a diagonal matrix representing transmission medium between transmitting unit 2 and said first one of said two antennas Λ_(t-2) ⁻¹ is a diagonal matrix representing transmission medium between said transmitting unit 1 and a second one of said two antennas Λ₂₋₂ is a diagonal matrix representing transmission medium between said transmitting unit 1 and said second one of said two antennas.
 15. A method carried out in a receiver for decoding received frame signals of a unit that transmits over p antennas, comprising the steps of: converting each received frame signal to frequency domain; in groups of p consecutive converted frame signals, combining said converted frame signals to form p intermediate signals; multiplying said intermediate signals by values related to transfer characteristics between said p antennas and said receiver, to obtain thereby equalized signals; converting said equalized signals to time domain, to obtain time domain estimate signals; and carrying out a decision regarding information symbols transmitted by said unit, based on said estimate signals.
 16. The method of claim 15 where said combining is linear combining.
 17. The method of claim 15 where said transfer characteristics employed in said step of multiplying are frequency domain characteristics of transmission channel between said p antennas and said receiver. 